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ABSTRACT 

We present high spatial resolution (5.1" x 2.7") Owens Valley Radio Obser- 
vatory (OVRO) millimeter array observations of HCN (J=l-0) and HCO + (J=l- 
0) emission in the inner 3 pc of the Galaxy (0.04 pc ~ 1"). The HCN and 
HCO + emission of the Circumnuclear Disk (CND) is distributed in a well-defined 
ring with a peak at a radius of 1.6 pc. The observed radial velocities are generally 
consistent with rotation at ~ 110 km s _1 (except along the western edge of the 
CND). The HCO + /HCN emission ratio is typically ~ 0.4 but with significant 
variations. The variations in the HCO + /HCN emission and absorption ratios 
can be attributed to greater abundances of HCO + in lower density regions both 
within the CND and along the line-of-sight. The HCN emission is well correlated 
with the H2 emission at 2.12 /mi both in the main emission lobes of the CND 
and also in four H 2 and HCN filaments. Multiple areas of interaction between 
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the ionized gas and the CND are also seen - the western arm of the minispiral 
is spatially and kinematically consistent with being the ionized inner edge of the 
CND, and the northern arm may be connected to the CND northeast extension. 

With the enhanced spatial resolution of the HCN maps, we resolve nu- 
merous dense molecular gas cores within the CND with characteristic diame- 
ter ~ 7" (0.25 pc). For 26 of the more isolated cores, we have measured sizes, 
velocity widths, and integrated fluxes. From these properties we estimated 3 
masses for each core: a virial mass assuming the cores are gravitationally bound, 
an optically thick mass from observed column densities of HCN, and a lower 
limit mass assuming the HCN emission is optically thin and shock-excitation is 
negligible. The virial and optically thick masses are in good agreement with a 
typical mass of (2-3) x 10 4 M Q and a total CND gas mass of 10 6 M . The internal 
densities implied by these core masses (assuming a uniform density distribution 
for each core) are on average (3-4) x 10 cm" 3 . The core densities are high enough 
to be stable against tidal disruption from Sgr A* and the central stellar concen- 
tration. This tidal stability suggests a longer lifetime for the CND. The high 
densities and masses within the cores might support star formation either in the 
CND itself or within a core infalling towards the inner parsec, thus providing 
a mechanism for the formation of the young stellar population observed in the 
inner arcseconds of the galaxy. 

Subject headings: Galaxyxenter — ISM :kinematics and dynamics — ISM:molecules 
— radio continuum:ISM — radio lines:ISM — stars:formation 



1. Introduction 

At a distance of only 8 kpc (Reid 1993), the Galactic Center presents a unique op- 
portunity to study in exquisite detail a supermassive black hole and its interactions with 
the surrounding medium. Extensive observations of the inner 3 parsecs of the Galaxy have 
revealed an extremely complicated, highly interacting region. Found within these inner par- 
secs are: a 4.0 x 10 6 M black hole (Ghez et al. 2005) coincident with the non-thermal 
radio source Sgr A* (Backer & Sramek 1999; Reid et al. 1999); a dense stellar population, 
including a central cluster of bright, young (less than a few Myr) emission-line stars (Morris 
& Serabyn 1996; Krabbe et al. 1995); X-ray emission, both in hundreds of point sources 
and in diffuse emission (Muno et al. 2003); streamers of ionized gas and dust, dubbed the 
"minispiral" (Ekers et al. 1983; Lo & Claussen 1983; Lacy et al. 1991; Scoville et al. 2003); 
and a clumped ring of molecular gas, the Circumnuclear Disk (CND), surrounding Sgr A*, 
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the minispiral, and the central star cluster (e.g. Becklin et al. 1982; Guesten et al. 1987; 
Wright et al. 2001). 

The CND was first detected by Becklin et al. (1982) as double- lobed emission at 50yum 
and 100/xm; it has subsequently been observed extensively at radio to infrared wavelengths 
(see Vollmer & Duschl (2001)). These observations indicate that the CND is a ring-like 
structure of molecular gas and dust with an inclination of 50° — 75° and a position angle 
of ~ 25° (Jackson et al. 1993). The ring is nearly complete in HCN but has noticeable gaps, 
especially to the north (Wright et al. 2001). The inner radius of the CND is well-defined at 
~ 1.5 pc (deprojected). Early single dish observations traced HCN, CO, and HCO + out to 
more than 7 pc (Genzel et al. 1985; Lugten et al. 1986; Serabyn et al. 1986; Guesten et al. 
1987), but a recent interferometric study suggests an outer edge at 3 - 4 pc (Wright et al. 
2001). The CND has a thickness of ~ 0.4 pc at the inner edge (Guesten et al. 1987; Jackson 
et al. 1993) and may expand to 2 pc at larger distances from Sgr A* (Vollmer & Duschl 
2001). CND kinematics are consistent with rotation at 110 km s^ 1 (Marr et al. 1993) with 
deviations in some regions, most notably to the west. These deviations are consistent with 
a warped disk created from the collision of two large clouds (Guesten et al. 1987) or a ring 
of multiple dynamically distinct streamers (Jackson et al. 1993). 

To understand the behavior of the CND in the complexity of the Galactic Center, and 
in particular to evaluate its stability, lifetime, and potential for star formation, it is critical 
to constrain the physical conditions of the CND molecular gas. Previous observations have 
found a total mass within the CND of only ~ 10 4 M (Genzel et al. 1985) and indicated 
densities of ~ 10 5 - 10 6 cm~ 3 (Marr et al. 1993). These densities are sufficiently low that 
the gas cores would be unstable to tidal disruption from Sgr A* and the central stellar 
population. However, Jackson et al. (1993) calculate densities within the CND cores of 10 6 - 
10 8 cm~ 3 from models of the HCN emission, noting that such densities would significantly 
increase the expected core lifetime. They do point out that modeling of other molecular 
and atomic tracers yields lower densities in keeping with the values of Marr et al. (1993). 
Here we present new observations of the CND in HCN (J=l-0) and HCO + (J=l-0) at a 
resolution of 5.1" x 2.7", a factor of 2-3 gain in spatial resolution over previous studies. 
With this enhanced resolution, we resolve substantial core structure within the CND and 
calculate typical core densities of 10 7 - 10 8 cm -3 , in agreement with the results of Jackson 
et al. (1993). In addition,, we measure typical gas core masses of a fewx 10 4 M , larger 
than previous estimates. 
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2. Observations 

The observations described here were taken from 1999 November to 2000 April at the 
Owens Valley Radio Observatory (OVRO) millimeter array. We obtained 10 tracks in four 
configurations (Table 1), resulting in a naturally-weighted synthesized beam FWHM of 5.1" x 
2.7". The observations used a 10 pointing mosaic to cover a 120" x 210" (4.8 x 8.4 pc) region 
centered on Sgr A*. Because of the low declination (-29°) of Sgr A*, each track was only 
3-4 hr in length. In total, 115 minutes were spent on each pointing during the observations, 
resulting in a la noise rms of ~ 30 mJy beam" 1 in each 6.7 km s _1 spectral channel and 
6 mJy beam -1 in the 1 GHz wide continuum band. 

We simultaneously imaged the HCN (J=l-0) and HCO + (J=l-0) transitions in the 
upper sideband at 88.63 and 89.19 GHz, respectively. HCN (J=l-0) and HCO+ (J=l-0) are 
both tracers of high density (> 10 5 -10 6 cm -3 ) molecular gas and are therefore ideal molecules 
for probing the dense CND cores. Each line was observed in 64 channels of 2 MHz (6.8 km 
s -1 ) width centered on Vlsr = km s -1 . The total velocity range was -213 to 206 km 
s -1 , easily covering the -150 to 150 km s _1 region where emission was detected in previous 
observations. 

In addition, we obtained a single track in the L configuration of the H 13 CN line at 86.34 
GHz. We observed the two pointings covering the northeastern and southwestern emission 
lobes of the CND and used these observations to estimate the HCN optical depth (thcn) for 
the CND (§4.2.2). 

A bright quasar was observed twice in each track for passband and flux calibration. The 
quasar fluxes, derived from comparisons with observations of Uranus and Neptune in nearby 
tracks, varied slowly over the six months of the observation by as much as 20% (Table 1), 
but this variation was accounted for in the individual flux calibration of each track. We 
estimate a total uncertainty in the flux measurements of ~ 10%. NRAO 530 was observed 
every 15 minutes for gain calibration, and based on the NRAO 530 observations, periods of 
low coherence on each baseline were removed. All of the calibration and editing was done 
with MMA (Scoville et al. 1993), while all cleaning and mapping used the Miriad software 
package (Sault et al. 1995). 

The dominant noise in these interferometer images is misplaced flux of bright emission 
peaks resulting from residual atmospheric and instrumental phase fluctuations. Using the 
strong signal from the continuum source and the closure phase relationship (see Cornwell 
& Wilkinson (1981)), the random phase fluctuations are corrected. Such iterative phase 
self-calibration was performed on the five pointings containing Sgr A* in the primary beam 
(diameter ~ 90" for 89 GHz) using line free channels to produce the continuum model by 
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cleaning down to the first negative clean component. The continuum model is dominated 
by point source emission from Sgr A*. The typical time interval for self-calibration was 30 
minutes, corresponding to the period between repetitions of the same map position. The self- 
calibration reduced the noise in the images by 20 - 30% and increased the total recovered flux 
for Sgr A*by ~ 50%. Since the self-calibration was based on line-free channels and contained 
only Sgr A*in the model, the CND emission line fluxes were increased negligibly by the self- 
calibration (typically by only 5% with a maximum increase of 15%). As noted above, the 
self-calibration was used to reduced the noise due to mis-calibrated and/or scattered flux 
of Sgr A*. Therefore, the lack of self-calibration on the five positions at the northern and 
southern edges of the field without Sgr A* in the primary beam did not significantly affect 
our determination of the CND line flux. 

The 10 pointings were mosaicked in the UV plane using the Miriad invert command 
with the "systemp" option, which gives low statistical weights to visibilities with high system 
temperatures. The maps were cleaned to a 1.5<r level in each individual channel using Miriad 
routine mossdi2, which is designed for mosaicked observations and cleans to a given level 
of the theoretical noise at each point in a mosaic rather than to a set flux level. We then 
subtracted the 85.91 GHz continuum map (Fig. 1) in the image plane from each channel 
map to produce the line emission maps. 

3. Results 

3.1. Sgr A* Continuum Observations 

The 85.9 GHz continuum map (Fig. 1) traces the ionized thermal emission of the min- 
ispiral and the non-thermal emission of Sgr A*. We adopt a location for Sgr A* at the peak 
of the 85.9 GHz continuum emission (17:42:29.3, -28:59:18.7, B1950). All the features of the 
minispiral are detected at > 3a level and are labeled in the figure. The peak continuum 
emission is 1.6 Jy beam -1 in the 5.1" x 2.7" beam. Wright et al. (2001) measure a peak flux 
of 2.14 Jy beam" 1 for a 13" x 4" beam. Convolving our observations to match their resolu- 
tion, we obtain a peak flux of 2.02 Jy beam" 1 , consistent with their results. For an aperture 
of 30" radius centered on Sgr A*, our map contains 5.38 Jy of continuum emission, compared 
with 8.9 Jy from Wright et al. (2001). This difference is likely due to the shorter baselines 
available in the Berkeley-Illinois-Maryland Association (BIMA) data, which enable detection 
of a greater fraction of the spatially extended components of the thermal emission. Tracing 
recovered flux as a function of distance from Sgr A*, we find very little additional flux recov- 
ered in apertures larger than 30-40". Nobeyama 45-m single dish observations (Tsuboi et al. 
1988), free of the missing short spacings problem, likewise observe little emission beyond a 
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30-40" radius from Sgr A* and measure 27 Jy within 60" of Sgr A* (Wright & Backer 1993), 
compared with ~ 5.4 Jy from our observations, suggesting we recover ~ 20% of the total 
3mm continuum flux. 



3.2. Line-of-sight Absorption 

Absorption by molecular gas along the line of sight to the Galactic Center can signif- 
icantly effect the observed HCN and HCO + emission at some velocities. Figure 2 shows 
HCN and HCO + spectra extracted from a 4" diameter region centered on Sgr A*. Both 
spectra show broad absorption from -80 to +70 km s _1 from Galactic center gas and nar- 
rower absorption features at -135, -50, and -30 km s _1 due to foreground Galactic arms at 
radii of 250 pc (Scoville 1972; Liszt & Burton 1978; Bieging et al. 1980), 3 kpc (Oort 1977), 
and 4 kpc (Menon & Ciotti 1970) respectively. 

Although the velocities of the HCN and HCO + absorption features are generally similar, 
the HCO + absorption is deeper than that of the HCN at -135 and 20 - 40 km s _1 . These 
deeper HCO + features suggest regions of lower density gas along the line of sight where the 
fractional abundance of HCO + has been enhanced because of deeper penetration by cosmic 
rays (Marr et al. 1993). 

3.3. HCN and HCO+ Emission from the CND 

The integrated HCN and HCO + intensity maps (Figs. 3 and 4) clearly show a ring-like 
CND with a gap in emission to the north of Sgr A* and reduced emission to the east. The 
overall CND gas distribution is similar to that seen in previous imaging, but these obser- 
vations detect more significant emission to the east of Sgr A* and because of the enhanced 
spatial resolution reveal much greater internal structure within the CND ring. The HCN flux 
is on average 2-3 times the HCO + flux (see §3.4). The CND flux distribution as a function 
of distance from Sgr A* (deprojected assuming an inclination angle of 60°) has a sharp inner 
edge at 40" (1.6 pc), an equally steep decline around 50" (2.0 pc), followed by a slowly declin- 
ing, non-zero flux out past 150" (Fig. 5). The sharp inner and outer edges strongly suggest 
the CND has a ring-like morphology (Wright et al. 2001). Despite its ring-like nature, we 
nevertheless maintain the traditional nomenclature of circumnuclear disk (CND). 

Channel maps of the HCN and HCO + emission binned in 27 km s" 1 (four channel) 
intervals are shown in Figures 6 and 7. The rotation of the CND is clearly evident, as 
the dominant negative velocity emission arises southwest of Sgr A* and loops around to the 
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northeast of Sgr A*, where the emission is primarily at positive velocities. The channel maps 
also clearly show locations with multiple velocity components. Perhaps the most prominent 
example is 20" west and 20" south of Sgr A*, where significant emission is seen around 
-40 and 70 - 100 km s" 1 . 

The rotation of the CND is additionally apparent in the flux-weighted average velocity 
map (Fig. 8). The velocity structure is continuous around the ring with the greatest neg- 
ative (approaching) velocity near the southwestern emission lobe and the greatest positive 
(receding) velocity near the northeastern emission lobe. The sharp velocity gradient along 
the western side of the CND is a result of the spatial coincidence of -50 and 70 - 110 km 
s -1 emission and the presence of strong absorption at velocities near km s _1 . The region 
to the east of Sgr A* has a velocity near km s" 1 and also suffers from large amounts of 
absorption, which explains the reduction in observed emission. We note that the primary 
velocities expected in the northern gap (~ 50 — 70 km s _1 ) do not suffer from severe line 
of sight absorption, thereby suggesting a true dearth of molecular emission in the northern 
gap. 

3.4. Comparing HCO+ to HCN Emission in the CND 

The HCO + /HCN abundance ratio is ~ 1 for typical Galactic molecular clouds (Nyman 
1983; Vogel & Welch 1983; Blake et al. 1987). Previous observations of the CND found that 
the HCO + emission traces the same regions as the HCN emission, but with HCO + /HCN as 
low as 0.06 - 0.2 (Marr et al. 1993). As seen in the integrated maps (Figs. 3 and 4), the 
relative strength of HCN versus HCO + emission in our observations is largely consistent with 
these previous measurements. For locations with > 3a emission in both HCN and HCO + , 
the mean HCO + /HCN ratio is 0.74. However, this result is misleading because we are 
only considering points with 3a emission in both tracers. There are almost 3 times as many 
locations with > 3a HCN emission without 3a HCO + emission as there are locations with 3a 
HCO + emission and < 3a HCN emission. Therefore, the true average HCO + /HCN emission 
ratio is significantly lower. Comparing the ratio of total HCO + to total HCN emission in 
each channel, we find a typical channel HCO + /HCN ratio of 0.4 (Fig. 9). 

Examining the variation of the HCO + /HCN ratio with velocity we find that the lowest 
ratios are found in the channels corresponding to the northeastern and southwestern emission 
lobes of the CND with the highest ratios around km s _1 (Fig. 9). This region of enhanced 
HCO + emission relative to the HCN emission corresponds roughly to the broad absorption 
trough around km s -1 (Fig. 2). While there are velocity ranges in which HCO + ab- 
sorption on Sgr A* is more significant than HCN absorption, especially around -140 and 
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20-40km s _1 (Fig. 2), the HCO+/HCN ratio in these regions is neither significantly enhanced 
nor reduced from the values in surrounding channels, suggesting that the differential line of 
sight absorption between HCN and HCO + does not significantly affect the analysis of the 
HCO + /HCN ratio in each channel. As we will see in §4, the southwestern and northeastern 
lobes include some of the most massive cores, while cores in our sample that are found at 
velocities around km s -1 tend to be the least massive. This suggests a trend where the 
lowest HCO + /HCN ratios correspond to high mass regions while the highest HCO+/HCN ra- 
tios are found in the lower mass regions. We estimate the mass of each core (see §4.2) in 
three independent ways. Two of these techniques do have a dependence upon the level of the 
HCN emission (or the abundance of HCN relative to H 2 ), which potentially biases us towards 
high mass cores in regions with high HCN flux (and thus low HCO + /HCN ratios). However, 
the third tracer, requiring an assumption that the cores are virialized, is independent of the 
observed HCN flux (and as we shall see provides density and mass estimates in good agree- 
ment with estimates from assuming the cores to be optically thick to HCN emission), thus 
strengthening the support for the correlation of low HCO + /HCN ratios in high mass regions. 
This correlation could result from the CND having enhanced HCN emission in more massive 
regions because of increased shock processing of the gas (perhaps in core-core collisions). 
Alternatively, if the lower mass regions have less shielding from cosmic ray ionization than 
do high mass regions, the fractional ionization rate would be higher in the low mass regions, 
leading to an enhanced HCO + abundance relative to HCN (Marr et al. 1993). 

One significant region of HCN emission without significant HCO + emission is the linear 
filament, the swath of HCN emission approximately 20" west and 35" north of Sgr A* (see 
location in Fig. 10). We discuss this feature in more detail below (§3.5) and argue that 
the HCN abundance may be enhanced relative to HCO + because of shocks in the linear 
filament. In contrast, there is a single region where HCO + emission dominates that of HCN. 
Between 17 and 30 km s" 1 in a region centered at 45" east and 5" south of Sgr A*, the 
median HCO + /HCN ratio is 1.22. The molecular gas cores found within this region tend 
to be some of the lowest mass and density cores in the CND, supporting the argument of 
higher HCO + abundance relative to HCN in lower mass and density regions due to enhanced 
cosmic ray ionization. 



3.5. CND and H 2 Emission 

In Figure 10 we overlay our integrated HCN emission map in contours upon a gray- 
scale image of the H 2 (1-0) S(l) emission at 2.12yum from Yusef-Zadeh et al. (2001). H 2 
can be excited by shocks or UV fluorescence. The H 2 emission is most prominent in the 
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CND in the northeastern and southwestern lobes (as a likely result of limb brightening) 
where the integrated HCN emission is also the greatest. Yusef-Zadeh et al. (2001) argue 
that the most likely excitation sources for the H 2 emission in the CND are UV radiation or 
the dissipation of energy from collisions of dense molecular gas cores within the CND. This 
core collision mechanism requires core densities of > 10 6 cm~ 3 , confirmed by our density 
measurements below (§4.3). Core collisions and/or shocks could also create regions with 
low HCO+/HCN ratios, as shocks can help fuel HCN emission and additionally increase 
HCN abundance relative to HCO+ (Mitchell & Deveau 1983). 

The linear filament is strongly detected in H 2 emission with two primary emission peaks 
coincident with two HCN cores (cores X and Y in Fig. 13 below). This feature is ex- 
tremely narrow in H 2 (~1") and contains an OH maser near the northern end. The emission 
may arise from shocks generated by the supernova shell of Sgr A East hitting the 50 km 
s _1 molecular cloud (Yusef-Zadeh et al. 2001). Such shocks may explain the extremely low 
HCO + /HCN ratio in the linear filament. 

A similar situation exists for the outer filament, which begins at the northern end of 
the CND northern arm. Weak HCN emission is seen coincident with the H 2 emission of the 
outer filament, but there is no accompanying HCO + emission. We note that this emission is 
weaker and does not appear in contours in the integrated emission HCN map in Figure 10 
but peaks around 10-15" W, 90" N. In addition, the outer filament coincides spatially and 
kinematically (peaking around 50 km s _1 ) with [Sill] emission at 34.8 fim (Stolovy 1997). 
The [Sill] emission and the HCN emission without accompanying HCO + flux both suggest 
a shock excitation mechanism for H 2 . 

Two additional regions of good spatial agreement between the H 2 and HCN are the 
northeastern and southern extensions. Both of these regions are unshielded from the radia- 
tion of the central stellar cavity, suggesting that the H 2 emission is caused by UV excitation. 
We discuss the northeastern extension in greater detail in Section 3.6. 

3.6. CND and the Minispiral 

To understand the relationship of the CND to the ionized minispiral, we utilized two 
maps of the ionized emission. The first, a Hubble Space Telescope (HST) NICMOS Paa 
map (Scoville et al. 2003) has an exquisite spatial resolution of 0.2". The second, a Very 
Large Array (VLA) H92a radio recombination line map (Roberts & Goss 1993) combines 
1" spatial resolution with 14 km s -1 spectral resolution to provide kinematic detail of the 
inner parsecs. 
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Figure 3 shows our integrated HCN emission map in contours on the color-scale Paa 
image. The most obvious interaction is along the minispiral's western arc, where the ionized 
gas is consistently immediately interior to the CND. The western arc has been suggested 
as the inner edge of the CND, ionized by the strong UV radiation field of the stars in the 
central cluster (Lo & Claussen 1983; Genzel et al. 1985; Serabyn & Lacy 1985; Guesten et al. 
1987). We confirm this relationship kinematically through examining the velocity centroids 
of the ionized and molecular emission along the western arc (Fig. 11). The velocities are 
in close agreement, particularly at a position angle of approximately -145°, at which both 
the HCN and H92a emission peak. While some overlap between the western arc and the 
CND exists, the ionized emission always lies interior to the peak HCN emission. The high 
column densities that we measure for the HCN emission along the western side of the CND 
correspond to visual extinctions of up to 700 magnitudes (Scoville et al. 2003), suggesting 
that in regions where the ionized emission observed in Paa overlaps the CND emission, the 
ionized emission must lie near the foreground side of the molecular gas. 

In addition, Figure 3 suggests that the peak emission from the northern arm of the 
minispiral may connect with the northeastern extension from the CND at ~ 5" E, 20" N. The 
enhanced spatial resolution of our observations reveals the northeast extension to be a longer 
and more sharply collimated feature than previously observed. The strong collimation in the 
molecular and ionized gas can be explained by pressure of mass-loss winds from emission-line 
stars in the inner parsec (Scoville et al. 2003). At the intersection between the northern arm 
and the northeastern extension the velocities are in remarkable agreement (Fig. 12). On the 
other hand, the northernmost parts of the northern arm appear to bend to the west slightly 
after meeting the northeastern extension and pass through the northern gap in the CND, and 
based on an increase in CND linewidths on either side of the northern gap and the presence 
of an OH maser within the gap (Yusef-Zadeh et al. 1999), Wright et al. (2001) suggest that 
the minispiral northern arm is infalling and creates the gap as it pushes through the CND. 
A definite conclusion on the relationship between the northern arm and the CND is not 
possible with these observations, but our improved spatial resolution allows us to suggest for 
the first time that the northeastern extension and the northern arm may be part of a single, 
highly collimated feature. 



4. Cores within the CND 

The enhanced spatial resolution of our HCN and HCO + observations resolves a great 
deal of internal gas structure within the CND, enabling us to analyze the molecular gas 
properties there using techniques similar to those used to derive the properties of Galactic 
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giant molecular clouds (GMCs) and their cloud cores. 

4.1. Core Identification and Measurements 

From the HCN maps we have identified 26 bright, isolated emission features, each a dense 
molecular gas core within the CND. These features, labeled A-Z, are marked in Figure 13. 
We required that each core included in this analysis be isolated in both position and velocity 
space from other emission peaks down to at least 50% of the peak emission for the core. 
Included in this sample are the vast majority of the brightest emission cores that met our 
requirements for isolation, as well as a number of lower emission level cores; however, we 
do not claim completeness to any particular flux level but rather intend the sample to 
be representative of the range of bright cores found within the CND. Some peaks in the 
integrated HCN emission map (Fig. 3), such as that 10" west of core F, are too close spatially 
and kinematically and thus are not included in our core sample. There are additional emission 
cores at larger distances from Sgr A* (e.g., 5" W, 80" N in Fig. 8), but we restricted our 
sample to cores contained in the main CND emission.. In Table 2 we detail the properties 
of the 26 cores, including location, size, velocity center and FWHM, and mass and density 
estimates. Each of these core properties is discussed in greater detail below. 

For each of the cores, position-velocity strip maps were generated in the north-south 
and east- west directions to use for measurements of core size and velocity width. In addition, 
the strip maps allowed us to reduce line of sight blending and separate multiple kinematic 
components (such as core pairs K/L and S/T). Examples of these position- velocity maps are 
given in Figure 14 for cores B, I, J, O, and V. In each of these plots we indicate the adopted 
half-power size and velocity FWHM estimates by the dashed lines. The measured apparent 
size FWHMs were deconvolved for the instrumental spatial resolution (beamsize) using the 
Gaussian approximation (i.e. sizet rue = size apparcnt — sizef nstrumcntal ). The mean core FWHM 
(corrected for the instrumental resolution) is 8.8 x 5.1" or 0.34 x 0.20 pc. Approximating 
the average diameter of each core by the geometric mean of its north-south and east-west 
FWHMs, we find that the mean core diameter is 6.7" (0.25 pc). 

In Table 2, we also include the projected and deprojected distance of each core from 
Sgr A*. To deproject, we assumed an inclination angle of 60° for the CND (Jackson et al. 
1993). The median and mean deprojected core distances are 1.67 pc and 1.80 pc respec- 
tively, and the distribution of the 26 core distances follows very closely the deprojected flux 
distribution found in Figure 5. 

HCN (J=l-0) has three hyperfine lines with relative intensities of 1:5:3 and relative 
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velocities of -7.07:0:4.8 km s -1 (Park et al. 1999). This velocity separation is on the order 
of the width of each channel (6.8 km s _1 ). We modeled the increase in measured velocity 
FWHM due to the hyperfine structure, as well as the instrumental velocity resolution (taken 
to be one channel width) and found it to be a minor effect for velocity FWHMs > 20 km 
s _1 but substantial for smaller widths. The simulations were used to correct the measured 
velocity widths of the CND cores. We find median and mean velocity FWHMs for the 
26 cores of 28.1 and 26.1 km s _1 respectively, corresponding to approximately 4 channels 
in our observations. It should be noted that the 1:5:3 intensity ratio of the three HCN 
(J=l-0) hyperfine lines was derived for optically thin emission and so may not be the exact 
ratio applicable in this situation, because the HCN emission is optically thick. However, 
comparison with linewidths measured from our HCO + observations for a smaller sample of 
cores suggests that the linewidths measured from the HCN emission and contained in Table 2 
are in agreement with the HCO + linewidths to within 15 - 20%. 



4.2. Core Masses 

We use three independent approaches to estimate or constrain the masses of the molec- 
ular gas cores in the CND: virial estimates assuming the cores are gravitationally bound, an 
optically thick estimate obtained by integrating the gas column density over the projected 
area of each core, and a lower limit obtained under the assumptions that the HCN in the 
cores is optically thin and the contribution from shocks to the HCN flux is minimal. Each of 
these approaches has its own limitations (eg. the cores may not be in virial equilibrium and 
the HCN to H 2 abundance ratio is very poorly known for the CND). In addition, the optically 
thin and optically thick estimates are affected by the uncertainty in the role that shock pro- 
cessing plays in producing the observed HCN emission. Shocks, particularly from core-core 
collisions, could substantially enhance the HCN abundance and hence the observed HCN 
flux per unit mass of gas. If the true HCN abundance is in fact enhanced, our estimates of 
the core masses and densities under both the optically thin and optically thick assumptions 
would be overestimated. The optically thin calculations therefore can be considered a lower 
limit to the core mass and density only under the assumption that the shock contribution 
to the HCN emission is negligible. 

Keeping in mind these assumptions and their possible effects, we find that the virial and 
optically thick methods yield substantially higher masses and densities for the CND than 
previously derived from lower resolution HCN imaging. 
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4-2.1. Virial Masses 

Virial mass estimates for the cores were computed from the measured sizes and velocity 
widths using 

3 CM 

^ = -^^M^ = 200R pc AVL s -! M (1) 

where R pc is the radius (pc) of the core (taken as the geometric mean of the deconvolved 
half widths in the north-south and east-west directions) and AV is the velocity line-width 
(FWHM) in km s -1 . We model the core as a uniform density sphere, introducing the factor 
of | in eq. 1. 



4-2.2. Molecular Column Density (Vhcn — 4j Masses 

The core masses can also be estimated by integrating the molecular gas column densities 
over the surface area of each core. The optical depth of the observed HCN emission can be 
estimated from the observed flux ratio of HCN and H 13 CN emission. Using the lowest 
resolution but highest signal-to-noise ratio data (L configuration tracks in Table 1) of the 
HCN line and the additional L configuration observation of the H 13 CN line we obtain a mean 
HCN (J=l-0) optical depth (t H cn) of 4 for the southwestern and northeastern CND emission 
lobes (assuming a C/ 13 C isotopic ratio of 30 for the Galactic center). This measured optical 
depth is in excellent agreement with previous estimates for the CND (Marr et al. 1993). 

The total column density of H2 can be derived from the HCN (J=l-0) optical depth 
(thcn) by: 



tvt thcn Sttz/ 3 AV (t HC n/4) AV km „-! [T X /(50K)] 2 

A1-0 fo (gi/g ) c 1 - e^ Xhcn/10 9 

where T x is the excitation temperature of the HCN rotational levels, Xhcn is the abundance 
of HCN relative to H 2 , Ai_ = 2.38 x 10" 5 sec^ 1 is the Einstein A coefficient for the HCN 
J=(l-0) transition and f = pp 2 -, for which B = 44.316 GHz is the rotational constant for the 
HCN (J=l-0) transition. The CND dust temperature is between 20 and 80 K (Becklin et al. 
1982; Mezger et al. 1989). Since at the very high densities encountered in the CND, the dust 
and gas should be in thermal equilibrium, we have scaled the HCN excitation temperature 
to 50 K, the midpoint of the dust temperature range. 

Integrating the H 2 column density (~ n H2 2 R) over the surface area (ttR 2 ) of the core 
and multiplying by a factor of 1.36 to account for the expected He mass fraction within the 
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core, one obtains the total mass of gas M co i, 



M col = 6.5 x 10 2 (N H2 /10 24 cm- 2 ) R 2 lpc M, 



■0 



(3) 



The major uncertainties in the column density mass estimate lie in the adopted abun- 
dance of HCN relative to H 2 and in the excitation temperature. Here we have scaled the 
equations and used a value of Xhcn — 10~ 9 . This value has been derived for a number of 
high density cores in Galactic GMCs (Blake et al. 1987) and also for the CND by Marshall 
et al. (1995) from analysis of multi-transition HCN observations. We should note that much 
higher values for X H cn such as 8 x 1CT 8 and 2 x 1CT 8 have been adopted by Marr et al. 
(1993) and Jackson et al. (1993), respectively. However, Marr et al. (1993) considered only 
T x = 150 — 450 K, compared with 50 K adopted here based on the CND dust emission. Pre- 
vious works have taken a wide range of values for T x from 20 - 450 K (see Table 2 of Jackson 
et al. (1993)), although most atomic and molecular studies have taken a value of T x > 100 
K. Since the derived column densities and masses scale as T 2 X^^, N , adoption of the higher 
HCN abundance and the higher T x by Marr et al. (1993) largely offset each other in their 
effect on the derived masses. Shock processing of the molecular gas would also increase the 
value of X H cn, leading to an overestimate of the core masses. 



Lastly, we obtain a lower limit to the core mass under the assumptions that shock- 
excitation is negligible and the HCN emission is optically thin and thus every spontaneous 
decay from the J = 1 rotational level yields an escaping line photon. The observed line flux 
then translates directly into a total rate of HCN photons (p p hotons) by: 



However assuming a rotational temperature and HCN abundance, we can also express the 
total rate of HCN photons by: 



4-2.3. Optically Thin Masses 



Pphotons 



LyAz/ _ SAV(f)47Td 2 



(4) 




M C oreXHCN 



fiAi_ 



(5) 



Pphotons 



m H2 



where ulh 2 is the mass of a hydrogen molecule, and by equating eqs. 4 and 5, calculate a 
minimum H 2 mass for each core. The optically thin assumption is clearly not true based on 
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our measurements of t H cn — 4 for parts of the CND, but it does allow for an estimate of 
the lower limit to the mass (assuming negligible shock contribution). After correcting for 
the mass fraction of He, we obtain 



M op , m .=28.6 ^ A . V X — ;' T 9 f ° K) M e (6) 

(ahcn/10 y ) 

where Sj y AVk m s- 1 is the total emission line flux of the core. We have assumed a distance 
to the CND of 8 kpc. If there were significant shock processing, then the actual masses and 
densities would be lower. 



4.3. Derived Masses and Densities 



In Table 2 we list the estimated masses (virial, t H cn = 4, and optically thin) for the 
26 CND cores identified in Figure 13. In Figure 15 we plot the distribution of these core 
masses. In these calculations, we assume T x = 50 K and Xhcn — 10~ 9 . The virial and 
thcn = 4 mass estimates are in good agreement, ranging from 2000 to 1.36 xlO 5 M Q for the 
selected cores with a median mass of 1.65xl0 4 and 2.43xl0 4 M Q , respectively. The lower 
limit, optically thin masses are typically a factor of 10 lower. The total mass contained 
within the 26 cores is 5.5 xlO 5 , 7.90 xlO 5 , and 6.79xl0 4 M assuming virial, t H cn = 4, and 
optically thin densities respectively. The 26 cores included in our sample cover a significant 
fraction of the dense HCN emission from the CND; therefore we estimate a total mass within 
the cores of the CND of 10 6 M Q . We note that Shukla et al. (2004) estimated virial masses 
for a sample of nine cores in the CND based on separate 7 x 3" OVRO HCO + imaging range 
from 3.0 x 10 3 to 4.5 x 10 4 M , similar to our virial mass estimates. The close similarity 
between the virial and thcn = 4 mass distributions (Fig. 15) supports the assumption that 
the cores are in virial equilibrium. 

The internal densities of the CND cores can be obtained from their mass estimates by: 



o ro ^6 (M/1000 M o ) 3 
n H2 = 3.58 x 10 6 v ' — cm" 3 (7) 

R-O.lpc 

where we again assume a uniform density sphere. The density estimates for each core under 
the three different mass estimates are contained in Table 2. 

We graphically show the derived core densities as a function of deprojected distance from 
Sgr A* in Figure 16. The core densities are significantly larger than the 10 5 - 10 6 cm -3 pre- 
viously estimated for the CND (Genzel et al. 1985; Marr et al. 1993). Indeed our mean virial 
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and thcn = 4 densities are 3.8 xlO 7 and 5.0 xlO 7 cm -3 , respectively. All of the optically 
thin density estimates are > 10 6 cm -3 . Our density estimates are, in fact, similar to those 
(10 6 — 10 8 cm~ 3 ) derived by Jackson et al. (1993) from single dish, multi-transition HCN 
measurements; their estimates are obtained from analysis of molecular excitation and hence 
are entirely independent of the virial and column density (thcn = 4) techniques used here. 
Similarly, Shukla et al. (2004) find virial densities of 4 x 10 6 to 3 x 10 7 cm~ 3 for their sample 
of nine CND cores. As we will discuss later, the high densities of our large core sample have 
important implications for core stability in the tidal field of Sgr A* and for the possible role 
of these cores in the formation of stars in the inner parsecs. 

A final independent check on our mass and density estimates of the CND is provided 
by millimeter and sub-millimeter dust emission measurements from the CND (Mezger et al. 
1989; Dent et al. 1993). Mezger et al. (1989) detected the CND dust at A = 1.3 mm. After 
subtracting substantial radio continuum from the ionized gas, they estimate a total mass 
of cold dust (at 20 K) and gas ~ 2 x 10 4 M Q in the northeastern and southwestern CND 
lobes. In 450 jum James Clerk Maxwell Telescope (JCMT) observations at 7" resolution, the 
flux of the CND is approximately 100 Jy in a distribution that closely follows that of the 
HCN observations of the CND (compare Fig. 3 in this work to Fig. 10a of Dent et al. (1993)). 
For a mass opacity coefficient of 24 gr cm -3 at 450 fim (Hildebrand 1983) this implies a mass 
of ~ 1 x 10 4 M for T d = 50 K. These estimates are more than a factor of 10 less than our 
measurements of the molecular gas mass within the CND. However, mass estimates based 
on observed dust emission require assumptions for the dust opacity coefficient and dust-to- 
gas abundance ratio - both of which could be quite different in the high density CND - 
and the assumption that the dust was optically thin at 450 /im. We note that the typical 
densities and sizes estimated in Table 2 for the cores yield typical mass column densities of 
50 to 100 gr cm -2 . For standard dust properties and abundances (Hildebrand 1983), these 
column densities imply optical depths of ~ 5 at 450/im. Thus the observed 450/im fluxes may 
considerably underestimate the CND mass because of opacity effects. Our 3 mm continuum 
maps show no evidence of dust emission associated with the CND molecular gas cores, but 
this is below our detection limit even for the very high column densities estimated above. 

One factor that could affect our measurements of the core masses and densities is the 
resolving out of substantial extended HCN flux. However, we compared our HCN maps to 
Nobeyama 45m HCN observations (Kaifu et al. 1987) and find that for the CND lobes we 
recover ~ 75% of the total HCN flux. 
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4.4. Tidal Stability 

The tidal shear on CND cores from the gravitational potential of Sgr A* and the central 
stellar population is extremely large. The total mass enclosed as a function of Galactic radius 
(L) can be modeled as: 

M G = 4.0 x 10 6 + 1.6 x 10 6 L^f M (8) 

This expression for total enclosed mass is taken from Vollmer & Duschl (2000) and modified 
to account for new measurements of the mass of Sgr A* (Ghez et al. 2005). 

Therefore the mean internal density required for a core to be tidally stable is : 

n H2 , ti dai = 2.87 x 10 7 (L- 3 + 0.4L pc L75 ) cm" 3 (9) 

We list the minimum densities for tidal stability (n H2i tidai) f° r eacn core i n Table 2. The 
mean value is 1.2 x 10 7 cm~ 3 . 

Comparing the calculated core densities with the minimum densities required for tidal 
stability, we find that the vast majority of cores are dense enough to be tidally stable. Indeed, 
all 26 cores have thcn = 4 densities sufficient for tidal stability, and 23 of the 26 cores 
have sufficient virial densities. We note that there is no correlation between measured core 
densities and deprojected distance from Sgr A* and that the virial and thcn = 4 densities 
are often 3-5 times greater than the minimal tidal stability density (Fig. 16). This is strong 
observational support for widespread tidal stability and virial equilibrium within the CND 
cores because if the cores were being tidally disrupted, their apparent virial densities would 
always approach the tidal stability limit as the cores gradually become tidally stripped. 

5. Discussion: Implications of CND Core Stability 

5.1. Lifetime of the CND 

The core densities that we measure for the CND are not only substantially higher than 
previous estimates but also on average 3-5 times larger than necessary for tidal stability. 
Thus the lifetime of the CND cores may be considerably longer than their internal dynamical 
timescales of ~ 10 5 yrs. 

The ionized gas within the minispiral has been estimated at 100 M Q (Sanders 1998). 
We calculate a total molecular mass within the CND cores of 10 6 M , giving an |Sg mass 
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ratio of ~ 10 4 . Estimates for the infall time of the ionized material are ~ 10 4 yr (Scoville 
et al. 2003; Vollmer & Duschl 2001; Sanders 1998). If even 10% of the molecular gas in 
the CND ultimately fell into the central parsec and became the ionized material now visible 
in the minispiral, the mass ratio between the CND and the minispiral and the minispiral 
infall time argue for a typical lifetime for the material in the CND of ~ 10 7 yrs. This time 
frame is a factor of 10-100 larger than previous estimates and is also much longer than the 
dynamical/orbital time for the CND. 

The neutral, atomic gas contained in the inner parsec has been estimated at 300 
M Q (Jackson et al. 1993). The bulk of this gas has been modeled as an infalling atomic 
cloud with the northern and eastern arms of the minispiral as the photoionized inner edges 
of this cloud. If this atomic mass is included in the mass conservation arguments above, 
this reduces the lifetime of the CND by a factor of 4, which is still significantly longer than 
previous estimates. 

The sharpness in the dropoff of the distribution of HCN flux and CND cores as a function 
of deprojected distance from Sgr A* (Fig. 5 and §4.1) has implications for the lifetime of the 
CND. More than 70% of the HCN emission from our observations is contained within the 
inner 2 pc. If the CND has a 10 7 yr lifetime, as suggested from the mass conservation 
argument, then there are sufficient molecular gas cores at radii beyond the CND that over 
10 7 yr could migrate to the CND, providing a means of maintaining its mass. Coil & Ho 
(1999, 2000) find evidence for such feeding of the CND from the 20 km s _1 molecular cloud 
at a 10 pc projected distance from Sgr A*and even detect cores within the streamer feeding 
the CND. In contrast, a short lifetime for the cores of the CND would suggest that the CND 
itself is a short-lived structure, as there is insufficient additional material (as traced either 
by the cores or by the raw HCN emission) outside the CND but within the inner few parsecs 
for replenishment. 



5.2. Formation of the Massive Young Stars Found in the Inner Parsec 

Several stellar populations are observed within the central parsec (Genzel 2001). These 
include a 1-10 Gyr stellar population dominated by emission from red giants (Genzel et al. 
2003), tens of blue supergiants tracing a very recent (2-7 Myr) burst of star formation (For- 
rest et al. 1987; Allen et al. 1990; Krabbe et al. 1991, 1995; Morris & Serabyn 1996; Tamblyn 
et al. 1996; Blum et al. 1996; Paumard et al. 2001), and a population of presumably young, 
dust-enshrouded stars (Becklin et al. 1978; Krabbe et al. 1995; Genzel et al. 1996, 2003). Gen- 
zel et al. (2003) describe the stellar population of the inner parsec as a metal-rich older star 
cluster, interspersed with a number of young stars, which dominate the emission in the inner 
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10". Scoville et al. (2003) detect a drop in the stellar surface brightness within 0.8" of Sgr A*, 
likely due to the interaction of the black hole with the surrounding stellar environment. 

Understanding the formation mechanism of the young stars within the inner parsec 
has proven extremely difficult. The enormous gravitational potential of Sgr A* prevents in 
situ formation without the aid of shocks or collisions unless molecular clouds have densities 
greater than 10 9 - 10 10 cm~ 3 . While we measure high densities within the cores of the CND, 
our measured densities are still 100-1000 times below this in situ formation threshold. On 
the other hand, star formation at greater distances (ie. > 10 pc) from Sgr A* and subsequent 
inward migration is also unfeasible because the time for two-body relaxation of a star to the 
inner arcseconds is greater than the observed age of the young stars. 

Our discovery of the high densities and masses within the CND cores provides a possible 
source for the young stars of the central parsec. While the transport of individual stars to 
the central region is impossible within the observed age of these young stars, star clusters 
can more quickly migrate to the inner parsec via dynamical friction. Portegies Zwart et al. 
(2003) estimate that a 6 x 10 4 M Q or greater cluster with a diameter of 0.3 pc or less could 
reach the inner parsec from a radius of 5 pc or less in the 2-7 Myr age of the young stars. 
This formation mechanism had been previously dismissed in part because of the lack of 
massive cloud cores that could produce such clusters in the inner 5 pc (Genzel et al. 2003). 
However, within our sample of 26 massive CND cores, almost all have diameters less than 0.3 
pes and masses greater than 10 4 M . Indeed the most massive core has a mass in excess of 
10 5 M . These most massive cores could potentially form a stellar cluster in the CND with 
high star formation efficiency and migrate to the Galactic center within 2-7 Myr, producing 
the observed young stellar population. For example, evidence suggests that IRS13 may be 
a remnant of one such cluster (Maillard et al. 2004). In addition, simulations of the CND 
suggest that the CND cores we observe are able to collapse and form stars before being 
shredded in the tidal field of the Galactic center (R. Coker 2004, private communication). 
Furthermore, numerical simulations have found that a 10 5 M Q cluster at the radius of the 
CND (or a less massive cluster harboring an intermediate mass black hole) could spiral into 
the inner parsec within the lifetime of the 2-7 Myr stars (Kim & Morris 2003; Kim et al. 
2004; Hansen & Milosavljevic 2003). Mid- infrared imaging of the CND (Telesco et al. 1996) 
has revealed significant emission at 10, 20, and 30 /zm, especially along the western arc of 
the minispiral and more diffusely over the entire western side of the CND. There are no 
significant peaks of mid-infrared emission coincident with the cores visible in the CND. Such 
a coincidence could have been an indication of active star formation. Sensitive high spatial 
resolution observations in the mid to far infrared (ie. 10-100/im) would allow for a more 
precise search for ongoing star-formation within the CND cores. 
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Cores within the CND reach densities of a few x 10 7 cm -3 , making them tidally stable 
at distances less than 1 pc from the Galactic center (Fig. 16). An alternate explanation, 
therefore, for the young nuclear star cluster is that a massive core remains tidally stable 
for the majority of its infall from the CND onto Sgr A* and forms stars during this infall. 
These new stars would then reach the inner parsec in less than 2-7 Myr. Depending on the 
efficiency of the star formation in the infalling core, multiple infall events might be necessary. 
These infall events could be spread over a period of a few Myr, explaining the observed age 
spread within the nuclear stellar cluster. 



6. Conclusions 

We have presented OVRO millimeter observations of HCN and HCO + emission from 
the CND at a spatial resolution of 5.1" x 2.7", a significant improvement in spatial resolution 
over existing observations. With this enhanced resolution, we were able to study in great 
detail the behavior of the CND in the gravitational potential of the Galactic Center and the 
possible role of the CND in the formation of the young stellar population of the inner parsec. 
In particular we find: 

1. The CND, as traced by HCN and HCO + emission is a well-defined structure with an 
inner radius of ~ 40" (1.6 pc) from Sgr A*and a sharp drop past ~ 50" (2.0 pc), suggesting 
a ring-like morphology. The main HCN and HCO + emission is consistent with rotation at 
~ 110 km s -1 , although emission at multiple velocities is found in some locations, particularly 
along the western edge of the CND. There is a lack of observed emission in the eastern part of 
the ring where radial velocities are near km s -1 because of strong line of sight absorption. 

2. The HCN flux is typically a factor of 2 - 3 greater than the HCO + flux within the 
CND. However, there is substantial variation in this HCO + /HCN flux ratio with density, 
likely because of a higher fractional ionization from cosmic rays and thus greater abundance 
of HCO + in lower density regions. In addition, shocks may play a role in enhancing the 
HCN abundance and flux. 

3. There is good spatial agreement between the H2 and HCN emission in the northeast- 
ern and southwestern CND lobes, as well as in the southern and northeastern extensions. In 
these regions, the H 2 is likely UV excited or excited by collisions of the dense CND cores. 
The linear and outer filaments are also seen significantly in H 2 and HCN (but not HCO + ) 
and are likely shock-excited. 

4. There are a number of interactions between the ionized gas of the minispiral and 
the CND. In particular, the western arc of the minispiral is both spatially and kinematically 
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consistent with being the ionized inner edge of the CND. In addition,, the minispiral northern 
arm may connect with the northeastern extension of the CND to form a single collimated 
structure. 

5. We identify a sample of 26 resolved molecular gas cores within the CND. These cores 
trace the overall HCN emission and have a characteristic diameter of ~ 6.7" (0.25 pc). The 
cores have a typical density of (3-4) xlO 7 cm~ 3 , and a typical mass of (2-3) xlO 4 M . The 
total mass within the CND is estimated at ~ 10 6 M . Some of the HCN emission may be 
enhanced by shock chemistry, implying that our estimates of the optically thick and optically 
thin core masses and densities (assuming a standard HCNabundance) are overestimates. 

6. The CND core densities are sufficient for tidal stability in almost all of the detected 
cores. Based on mass conservation arguments, a 10 7 yr lifetime for the CND is predicted. 

7. The core density estimates under the assumption of virialization and under the 
assumption of an optical depth thcn = 4 are similar and often 3-5 times the minimum density 
for tidal stability, suggesting that many of the cores are indeed gravitationally bound. 

8. The enhanced core densities and masses may explain the formation of the massive 
young (2-7 Myr) stars found in the inner parsec. The cores are tidally stable and can undergo 
star cluster formation while within the CND. Alternatively, these cores are stable even at 
closer distances to Sgr A* and could form stars during an infall process. In either scenario 
stars can be transported to the inner arcseconds within the observed age of the stars. 
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Table 1. Summary of Individual Observations 



Configuration 


Largest 


Date 


Usable 


Source 


3c273 


3c354.4 




Baseline 




Baselines 


Integration Time 


Flux 


Flux 




(m) 






(minutes) 


(Jy) 


(Jy) 


L 


115.0 


1999 Nov 13 


15 


114 




7.0 






2000 Apr 25 


14 


140 


9.0 


6.0 


E 


119.3 


2000 Mar 16 


15 


140 


9.1 


7.2 






2000 Mar 17 


15 


136 


9.1 


7.2 


H 


241.7 


1999 Dec 17 


15 


140 


9.8 


8.5 






2000 Jan 4 


8 


70 


9.8 


8.5 






2000 Jan 5 


15 


139 


9.8 


8.5 






2000 Jan 8 


15 


140 


9.8 


8.5 


U 


483.3 


2000 Feb 6 


15 


130 


9.5 


9.0 






2000 Feb 8 


15 


120 




9.0 



Table 2. Measured Core Properties 



(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) (14) 


(15) 


Core 


RA Dec 


Proj. 


Dcproj . 


Size 


Cent 


Vcl. 


Virial 


T HCN — 


4 Opt. Thin 


Virial 


thcn — 4 Opt. Thin 


Min Stable 


ID 


Offset Offset 


Dist. 


Dist. 


FWHM 


Vcl. 


FWHM 


Mass 


Mass 


Mass 


Dens 


Dens Dens 


Dens 




( from Sgr A*) 


(pc) 


(PC) 


(PC) 




(km/s) 




(X 10 3 


M Q ) 




(X 10 6 cm- 3 ) 





A 


9.2 


32.0 


1 


29 


1 


34 





17 


107 


30 


9 


16.55 


15.04 


0.90 


91.60 


83 


23 


4 


98 


18.83 


B 


10.8 


40.0 


1 


61 


1 


68 





27 


139 


31 


6 


27.33 


38.53 


1.85 


38.04 


53 


64 


2 


57 


10.70 


C 


24.0 


40.8 


1 


84 


1 


86 





34 


105 


43 


7 


65.63 


83.82 


8 


24 


46.37 


59 


22 


5 


82 


8.30 


D 


27.6 


34.8 


1 


72 


1 


86 





43 


101 


45 


5 


88.84 


135.96 


10.92 


32.29 


49 


41 


3 


97 


8.34 


E 


25.2 


26.8 


1 


43 


1 


62 





22 


78 


30 


8 


20.92 


24.32 


2 


85 


55.97 


65 


06 


7 


62 


11.59 


F 


48.8 


26.0 


2 


14 


3 


10 





28 


72 


20 





11.22 


25.49 


2 


55 


14.65 


33 


28 


3 


33 


2.55 


G 


50.0 


8.8 


1 


97 


3 


42 





29 


65 


13 


1 


4.98 


18.05 


1 


80 


5.76 


20 


87 


2 


08 


2.05 


H 


22.0 


-0.8 





85 


1 


61 





24 


-17 


17 


3 


7.10 


15.78 





74 


15.33 


34 


05 


1 


60 


11.87 


I 


22.0 


-10.0 





94 


1 


87 





26 


-18 


15 





5.77 


16.12 


1 


23 


9.68 


27 


06 


2 


06 


8.15 


J 


4.8 


-25.2 





99 


1 


42 





25 


-37 


11 


7 


3.41 


11.72 


1 


69 


6.41 


22 


01 


3 


16 


16.29 


K 


2.4 


-32.8 


1 


28 


1 


67 





25 


-71 


24 


5 


15.27 


25.67 


2 


06 


26.76 


44 


99 


3 


60 


10.81 


L 


1.6 


-33.2 


1 


29 


1 


66 





19 


-38 


12 


6 


3.06 


7.51 





78 


12.52 


30 


77 


3 


21 


11.02 


M 


-3.2 


-37.6 


1 


46 


1 


70 





26 


-64 


19 





9.49 


21.47 


2 


29 


14.78 


33 


43 


3 


56 


10.28 


N 


-16.4 


-43.6 


1 


81 


1 


82 





22 


-64.5 


14 


7 


4.78 


11.59 





93 


12.83 


31 


15 


2 


50 


8.73 


O 


-21.2 


-34.4 


1 


57 


1 


60 





33 


-108 


36 


5 


43.96 


64.58 


9 


23 


35.04 


51 


48 


7 


36 


12.04 


P 


-21.2 


-24.8 


1 


27 


1 


39 





21 


-73 


28 


2 


16.80 


20.49 


2 


34 


50.82 


62 


00 


7 


06 


16.93 


Q 


-26.8 


-20.0 


1.3 


1 


68 





37 


-38 


12 


6 


5.97 


28.69 


2 


72 


3.28 


15 


74 


1 


49 


10.69 


R. 


-18.0 


-7.6 


0.76 


1 


16 





14 


-37 


12 


2 


2.02 


3.66 





24 


22.91 


41 


63 


2 


77 


27.06 


S 


-23.2 


-6.0 


0.93 


1 


55 





32 


89 


34 


4 


37.56 


56.18 


2 


67 


33.81 


50 


57 


2 


40 


12.96 


T 


-24.0 


-5.6 


0.96 


1 


62 





24 


44 


33 


5 


26.39 


30.08 


1 


45 


58.23 


66 


36 


3 


19 


11.73 


U 


-16.8 


-5.6 


0.69 


1 


10 





24 


79 


34 


7 


29.01 


32.65 


1 


91 


59.71 


67 


20 


3 


93 


30.90 


V 


-10.8 


10.8 


0.59 


1 


13 





19 


58 


24 


2 


11.27 


14.59 


2 


42 


45.11 


58 


41 


9 


70 


28.95 


w 


-7.6 


23.2 


0.95 


1 


47 





22 


56 


27 


9 


17.25 


22.30 


2 


07 


45.28 


58 


52 


5 


44 


14.89 


X 


-24.4 


27.2 


1 


42 


2 


67 





16 


64 


28 


4 


13.04 


12.05 





75 


88.59 


81 


85 




12 


3.57 


Y 


-18.8 


39.2 


1 


69 









21 




35 


3 


25.90 







94 


91.92 




71 




97 


3.15 


Z 


-3.6 


40.4 


1 




2 


08 





24 


58 


39 


6 


37.63 


37.00 


2 




78.22 


76 


91 


4 


93 


6.33 


Mean Core 






1 


32 


1.80 





25 




26 


1 


21.20 


30.70 


2 


61 


37.92 


49.91 


4 


09 


12.26 



to 
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RA (B1950) 

Fig. 1. — Continuum emission at 85.9 GHz from Sgr A* and the minispiral. The contours 
start at 15 mJy beam" 1 and increase by factors of 2 2 . The main components of the minispiral 
are detected in the continuum image and labeled above. The beam FWHM is indicated by 
the filled oval. 
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Fig. 2. — HCN (solid) and HCO + (dashed) spectra extracted from a 4" diameter region 
centered on Sgr A*. Absorption features at -135, -50, and -30 km s -1 , as well as an absorption 
trough centered around km s -1 , are detected with varying depths in HCN and HCO + , 
perhaps indicating differing densities of absorbing material. 
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Fig. 3. — Integrated HCN emission map in contours overlaid upon the color-scale HST 
NICMOS image of extinction corrected Paa emission in the Galactic center region (Scoville 
et al. 2003). Coordinates are offset from Sgr A*. Emission contours (solid) are at [0.2, 0.4, 
0.6, 0.8, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0]*6.765 Jy beam" 1 km s" 1 while the 
absorption contours (dashed) are at [-0.5, -1.0, -2.5, -5.0, -10.0] *6. 765 Jy beam -1 km s _1 . 
The color-scale units, indicated on the right, are 10~ 16 ergs cm -2 s _1 per pixel. 
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Fig. 4. — Integrated HCO + emission map in contours overlaid upon the color-scale HST 

NICMOS image of extinction-corrected Paa emission in the Galactic Center region (Scoville 

et al. 2003). Coordinates are offset from Sgr A*. Contour levels are the same as for the 

integrated HCN emission map (Fig. 3). The color-scale units, indicated on the right, are 
1Q-16 er g g cm -2 s -i p er pj xe i 
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Distance frem Sqr A* (") 



Fig. 5. — Average integrated HCN flux per pixel as a function of distance from Sgr A*. 
Distances are marked in both arcseconds (bottom) and parsecs (top). The dashed line 
corresponds to the projected flux distribution, while the solid line is the deprojected flux 
distribution assuming an inclination angle of 60° (Jackson et al. 1993). Note that the area 
under each of the graphs is not identical because different numbers of pixels fell in the various 
distance intervals for the projected and deprojected distributions. 
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Fig. 6. — HCN emission binned into 27 km s" 1 (four channel) intervals. The lowest contour 
is at 5c (lei = 15 mJy beam -1 for four channels) and each subsequent level increases by 5a. 
The central velocity in km s -1 is listed in the top left corner of each map. 
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Fig. 7. — HC0 + emission binned into 27 km s _1 (four channel) intervals. The contour levels 
are identical to the levels in the HCN channel maps (Fig. 6), so the significantly reduced 
HC0 + emission relative to HCN is readily apparent. The central velocity in km s _1 is listed 
in the top left corner of each map. 
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RA offset (arcsec) 

Fig. 8. — Flux-weighted average velocity (moment 1) map of HCN emission in color with 
integrated HCN emission in contours. Only individual channel emission at > 2cr is included. 
The velocities are with respect to Vlsr = km s -1 , and negative velocities (shown in blue) 
correspond to approaching emission. Emission contours begin at 3.4Jy beam -1 km s _1 and 
are spaced by 6.8Jy beam -1 km s" 1 . 



-100 100 
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Fig. 9. — Ratio of total HCO + to total HCN emission in each channel as a function of 
velocity. Only 3a positive emission was included in the flux ratio measurements. The noise 
level is virtually identical for both the HCN and HCO + maps. 
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Fig. 10. — HCN emission in contours overlaid upon a reverse gray-scale H2 (1-0) S(l) line 
emission map (Yusef-Zadeh et al. 2001). Significant features both in the H 2 emission and 
the HCN CND emission are labeled, and Sgr A* is marked by a circle. The contours begin 
at 5.1 Jy beam" 1 km s" 1 and are spaced by 10.1 Jy beam -1 km s _1 . 
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Fig. 11. — Velocity centroids of HCN emission (line) and R92a emission (diamonds) along 
the western function of angle (degrees east of north) from Sgr A*. The region of peak 

emission for both the ionized western arc and the western portion of CND, corresponding to 
Core O in Figure 13, is indicated. 
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Dist N of Sgr A* (") 

Fig. 12. — HCN (solid line) and H92a (dashed line) scaled flux distributions along the 
minispiral northern arm and CND northeastern extension. Velocity centroids are shown for 
the HCN emission (+'s) and H92a emission (diamonds). The region of intersection between 
the northern arm and the northeastern extension is labeled. 
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Fig. 13. — HCN integrated emission contour map with locations of the isolated cores exam- 
ined in Table 2. Cores K/L and S/T are pairs of cores at approximately the same spatial 
location but distinct in velocity space. The HCN contours are the same as in Figure 3. 
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Fig. 14. — Sample of position-velocity cuts in the horizontal (left) and vertical (right) di- 
rections centered on five cores identified in Figure 13. The horizontal and vertical cuts are 
summed over 0.36" and 0.2" respectively. The contours are set at 12.5% increments of the 
peak intensity with the lowest green contour at 50% and the lowest red contour at 75%. The 
asterisk indicates the location of peak intensity, and the dashed lines mark the region used 
for estimating the velocity and size FWHMs for each core. 
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Fig. 15. — Distribution of virial (solid), t H cn = 4 (dot-dashed), and optically thin (dashed) 
masses in logarithmic mass bins. 
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Fig. 16. — Distribution of virial (asterisks), t H cn = 4 (triangles), and optically thin (dia- 
monds) densities as a function of deprojected distance from Sgr A*. The minimum density 
for tidal stability is shown by the solid curve. 



